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Abstract The magnetic susceptibility of the optimally 
doped Pri_ x LaCe 2: Cu04 in the superconducting state 
is calculated using the t-J model of Cu-0 planes, Mori's 
projection operator technique and the dispersion of elec- 
tron bands derived from photoemission experiments. 
The electron band folding across the antiferromagnctic 
Brillouin zone border, which is inherent in the crystal, 
leads to a commensurate low- frequency response. The 
same band folding causes the appearance of a supple- 
mentary spin-excitation branch. The coexistence of the 
two spin-excitation branches explains two maxima ob- 
served in the frequency dependence of the susceptibility. 
The calculated momentum and frequency dependencies 
are close to experimental observations. Similarities and 
differences in the magnetic responses of electron- and 
hole-doped cuprates are discussed. 
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1 Introduction 

Magnetic responses of n- and p-type cuprate perovskites 
are essentially different. In the former crystals, the low- 
frequency susceptibility is commensurate, while in the 
latter it is incommensurate, heaving peaks at momenta, 
which differ from the antifcrromagnetic (AF) wave vec- 
tor Q = (tt, 7r) [HI2] [hereafter, I use the two-dimensional 
(2D) notations for wave vectors in a Cu-0 plane with 
the lattice spacing set as the unit of length] . In n-type 
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cuprates the dispersion of the susceptibility maxima re- 
sembles a cone with the apex point at the frequency 
u> — and at the momentum Q [3]H3 . In p-type cuprates 
this dispersion has the hourglass shape with the waist 
at Q and at the frequency 25 — 50 meV in moderately 
doped crystals [2]. 

The susceptibility of n-type cuprates was calculated 
in the random phase approximation in Refs. [5,6 . The 
obtained result depends largely on the momentum de- 
pendence of the electron four-point vertex C/ q , which 
choice is, to a certain extent, arbitrary in this approach. 
It is known that C/q peaked at Q tends to suppress in- 
commensurability. Such a vertex was chosen in Refs. [SJ 

However, even with this vertex and with simple 
tight-binding approximations used for the electron dis- 
persion in these works the low-frequency susceptibility 
demonstrates weak incommensurability [5] . In the con- 
sideration of hole-doped cuprates the vertex is usually 
set to a constant (see, e.g., [7] and references therein). 
The application of such a vertex to n-type cuprates 
leads to pronounced incommensurability, which is in- 
consistent with experimental observations [5]. 

In this article, the t-J model of Cu-0 planes is 
used to calculate the dynamic magnetic susceptibility 
of Pri_ a; LaCe a ;Cu04 (PLCCO) in the superconduct- 
ing state. For x = 0.11 — 0.12 the magnetic response 
of this crystal was studied in a series of neutron scat- 
tering experiments [TPPPP[TuTin] . the electron dis- 
persion was derived |12j from photoemission data 13 , 
and estimates for the value of the superconducting gap 
were obtained from high-resolution scanning tunneling 
microscopy measurements [TTlll4| . The use of the t-J 
model allows one to take proper account of strong elec- 
tron correlations inherent in n- types cuprates. To calcu- 
late Green's functions constructed from Hubbard oper- 
ators of the model Mori's projection operator technique 
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[To] is used. For p-type cuprates this approach allowed 
us to reproduce the observed momentum and frequency 
dependencies of the susceptibility using hole dispersions 
derived from photoemission [16U17] . The cases of super- 
conducting and pseudogap phases were considered. 

The electronic structure of n-type cuprates is char- 
acterized by band folding across the AF Brillouin zone 
border [TUT2] . This fact plays a central role in the forma- 
tion of the commensurate low-frequency response. The 
momentum dependence of the low-frequency suscepti- 
bility is governed by the spin-excitation damping, which 
peaks at Q due to the band folding. This peculiarity of 
the electron-doped cuprates leads also to the appear- 
ance of a supplementary spin-excitation branch. As a 
consequence the dispersion of the susceptibility max- 
ima has the shape of a cone. Its upper part is formed 
by the nested into each other branches of the usual 
and supplementary spin excitations. The part near the 
apex point is determined by the spin-excitation damp- 
ing. The same two parts can be singled out in the dis- 
persion of hole-doped cuprates, with the difference that 
the supplementary spin excitations are lacking there 
and the low-frequency spin-excitation damping peaks 
at incommensurate momenta. In the frequency depen- 
dence of the susceptibility, the usual and supplementary 
spin excitations manifest themselves as two maxima or 
a maximum and a shoulder. Such frequency dependen- 
cies were recently observed in PLCCO [TUlIll) . 

2 Main formulas 

Formulas for the magnetic susceptibility of the 2D t-J 
model were derived in Ref. j!6j using Mori's projection 
operator technique. In this approach, the formally exact 
expression for the susceptibility can be obtained [15 
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Here k is the 2D wave vector, other parameters are 
expressed through correlators of the spin-i operators 
s k' s k an d their time derivatives. In the case of the t-J 
model these parameters read 



ftk = 4 + J|Ci|)(l-7k) 



(2) 



are correlators of the electron a k ^ and spin operators 
on the neighboring sites, a is the spin projection, the 
angular brackets denote the statistical averaging and N 
is the number of sites. The parameter a serves for cor- 
recting the decoupling procedure [TSlfTD] used for deriv- 
ing Eqs. © and (0). For small electron concentrations 
a ~ 1.7 [2D]. The quantities Wk and TT(kcj) are the fre- 
quency of spin excitations and their polarization opera- 
tor, which contains contributions from interactions with 
electrons and with other spin excitations. The param- 
eter 5 in Eq. ([3]) describes a gap in the spin-excitation 
spectrum near Q due to temperature fluctuations [T9] 
and/or the interaction with electrons [2DJ. The magni- 
tude of this gap is directly connected with the correla- 
tion length of the short-range AF order, and it grows 
with the electron concentration x w x and temperature 
T. In the case of hole doping the magnitude of the gap 
determines the frequency co r of the waist in the hour- 
glass dispersion of the susceptibility maxima [16] . 

Calculations are essentially simplified for zero tem- 
perature. For this case, the polarization operator reads 

n{\tu) = il x (kw) + n 2 {kuj) + 7T 3 (kw), (4) 
771 (M = lA^T E E [ 2 A 2 (kqq') + /I(kqq')] 

qq' tt' 
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16 J2a ^ 1 + 7$T\) (1 - 7k) {S + 1 + 7k) ' n *™ = mz EE /Kkq) 



(3) 

where to and J are the hopping and exchange constants 
between neighboring sites in the t-J Hamiltonian, 7k = 
[cos(fc x ) + cos(fcy)] /2, 

= T? E ^ k ( a L a ^) and c i = 4 E ^ ( s i s k) 
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iq and J q are Fourier transforms of the hopping and 
exchange constants in the t-J Hamiltonian. 

Equations (J!])-© were derived using the following 
normal AhQsoj) = — it" 1 Im((ak a -|aj c )) and anomalous 
j4i 2 (ku;) — 7r _1 Im((ak^|akj,)) electron spectral functions: 
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e± = i(e k ±e k - Q ), £ k =yJ(e^ + A}, 
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The dispersions i-Ekrj Eq. (fTU]) . correspond to elec- 
tron bands of a crystal with the superconducting gap 
function A k and with the band folding across the AF 
Brillouin zone border, which is characterized by the po- 
tential A f [T2ll2Tj . This potential and parameters of the 
initial tight-binding dispersion 

£k = — 2t [cos(k x ) + cos(fc y )] — At i cos(k x ) cos(k y ) 
— 2i 2 [cos(2fc x ) + cos^fey)] — 4i 3 [cos(2fca;) cos(fc y ) 
+ cos(fca;) cos(2fc l/ )] — 4^4 cos(2fc x ) cos(2fcj,) — /i (11) 

were obtained |12j by fitting the photoemission data 
Q3] in Pr .89LaCeo.iiCu0 4 : t = 0.12 eV, h = -0.06 eV, 
t 2 = 0.034 eV, t 3 = 0.007 eV, t 4 = 0.02 eV, fj, = 
-0.082 eV and A f = 0.14 eV. In the below calcula- 
tions, it was expected that only the nearest-neighbor 
exchange constant is nonzero, and its value J = 0.12 eV 

Equations ©-(0 describe four processes of the spin 
excitation transformation - its conversion to an electron- 
hole pair with assistance of another spin excitation [Eq. 
J5J], the decay into three spin excitations [Eq. (0] and 
the conversion into an electron-hole pair, direct and as- 
sisted by an electron [both processes are combined into 
Eq. |(7J|]. In the considered range of parameters this lat- 
ter term makes the main contribution to 77(kw). 

Notice that to does not coincide with the parame- 
ter f of the t-J Hamiltonian in Eqs. (0 and (0, since 
the interaction between electrons and spin excitations 
leads to a considerable reduction of the electron band- 
width. Based on the exact diagonalization results [35] 
I set to /to = 2 and supposed the same relation be- 
tween other unrenormalized and renormalized hopping 
constants in Eq. ([5]). Only the nearest-neighbor hoping 
constant was left in Eqs. and (0, since other kinetic 
terms introduce merely small corrections. Also a small 
third harmonic in the superconducting gap function [T] 
influences only slightly the susceptibility and the major 
part of calculations was carried out with the d-wave gap 
function A k = A[cos(k x ) — cos(k y )]/2. The gap param- 
eter A varied in the range 4.8 — 9.6 meV. These values 
are close to the estimates obtained by high-resolution 
scanning tunneling microscopy |11II14] . 

In Eq. (0, the parameter Z takes into account the 
difference of the operators a\^ a and a ka , which are the 
Hubbard operators, from the fermion creation and an- 
nihilation operators and the fact that in the case of 
strong electron correlations a sizable part of the spec- 
tral weight is concentrated in an incoherent part of the 
electron spectral function. The spectral weights of the 
coherent Z and incoherent Z' parts satisfy the relation 



Z + Z' = 
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Since the incoherent continuum does not produce peaks 
in the susceptibility, only the coherent parts of the spec- 
tral functions were taken into account in (J9j) . Supposing 
that in the n-type case the ratio ZjZ' is approximately 
the same as in the p-type case, Z was set to 1/6 for 
to = 2J. In the calculations, the infinitesimally small 
quantity 77 in Eqs. ([5])- ([7]) was substituted by the ar- 
tificial broadening r = 0.3 — 7 meV. Such values are 
usually used in calculations of the susceptibility. The 
parameter r influences mainly widths of maxima in 
the frequency dependence of the imaginary part of the 
susceptibility x"(kw). 

As mentioned above, the parameter 6 defines the lo- 
cation of the bottom of the spin-excitation dispersion 
iv r , which corresponds usually to the maximum in the 
frequency dependence of x"(Qcl>). The value of 6 was 
fitted to adjust the calculated maximum to the exper- 
imental one. Obtained values of S varied in the range 

2 x 1CP 4 — 4x 10~ 3 , which corresponds to u r of the order 
of several meV. These S and Lo r are much smaller than 
the respective parameters in moderately doped p-type 
cuprates, which is the consequence of the fact that the 
considered crystal is close to the boundary of the long- 
range AF order. In this connection, it is worth noting 
that in some works a weak static AF order was observed 
in the optimally doped PLCCO at low temperatures 
(see, e.g., Refs. [2111]). At the same time it was pointed 
out that low-frequency maxima are broader than the in- 
strumental resolution, and a finite correlation length of 
the AF order was obtained [3] , which implies a nonzero 
value of 5. The question about possible coexistence of 
superconductivity and long-range AF order is still open 
due to the influence of sample inhomogeneity . Fol- 
lowing results of Refs. [9,23] in this work it is supposed 
that the long-range AF order is lacking in the super- 
conducting phase and S is finite. 

3 Results and discussion 
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Fig. 1 The HWHM of the susceptibility maxima in reciprocal 
lattice units in PLCCO from Ref. g] (as = 0.11, T c = 25.5 K, 
T = 6 K, red squares), Ref. [3] (x = 0.12, T c = 21 K, T = 
7 K, blue triangles), and in the present calculations with the 
parameters given in the text and S = 0.001, A = 9.6 meV 
(black circles). Above u Rd 0.1 eV results of Ref. [3] give the 
dispersion of resolved incommensurate peaks. 
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Fig. 2 The momentum dependence of x"(kw) for u> = 
4.8 meV (a), 16.8 meV (b), 36 meV (c), 48 meV (d) and 
108 meV (e). The wave vector varies along a diagonal of the 
Brillouin zone, k = (k,k). Parameters are the same as in 

Fig.m 



The dispersion of the susceptibility x"(kw) at half max- 
imum, obtained in the present calculations, is shown 
in Fig. [TJ This quantity can be compared with results 
of Refs. [2,3,4 ], in which an unresolved commensurate 
peak is observed in momentum cuts up to u ps 0.07 eV 
[3] or even up to 0.17 eV 011] • Therefore, in these 
works frequency dependencies of the half-width at half- 
maximum (HWHM) of these peaks are given. These 
dependencies are also shown in Fig. \T\ As can be seen 
from this comparison, the calculations reproduce the 
experimentally observed cone-shaped dispersion with 
the apex point at lu — and k=Q. Moreover, the ob- 
tained q widths of peaks are in reasonable agreement 
with those observed in Refs. [HH]. This agreement is 



not a result of a r fitting - the q widths depend only 
weakly on this parameter. As observed in Ref. [3], above 
to w 0.1 eV the commensurate maximum splits into in- 
commensurate peaks, which disperse like spin waves. 

The calculated momentum cuts of x"(kw) along a 
diagonal of the Brillouin zone are shown in Fig. [5] For 
small frequencies the susceptibility peaks at Q [Fig. 
[21(a)] - the magnetic response is commensurate. This re- 
sult is retained for other used values of A and 8. As fol- 
lows from the above formulas, the main reason for such 
a response is the electron band folding across the AF 
Brillouin zone border. This folding produces nested low- 
frequency equi-energy contours with the nesting vector 
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Fig. 3 The equi-energy contours of electron dispersion (|10|l . 
(ITU for £j = ±6 meV and A = 3.6 meV (black ovals). The 
blue dashed box is the boundary of the Brillouin zone. Red 
arrows show transitions making the main contributions to 
the polarization operator Q for 8 = 0.001, uj = 1.2 meV and 
k = Q. 

Q (see Fig. [3]). Transitions between electron states on 
these contours make the main contribution to polar- 
ization operator (0. As a consequence — ImTJ(kw) in 
the numerator of the formula for x"(ku;) peaks sharply 
at Q. With a weaker momentum dependence of the 
denominator at small oj, the numerator controls the 
behavior of the susceptibility. The denominator also 
contributes to the appearance of the commensurate re- 
sponse, reaching the minimum at Q. Notice that in p- 
type cuprates the situation is different - in this case 
— Im77(kw) has sharp peaks at incommensurate mo- 
menta, ft is these peaks that produce the down-directed 
branch of the hourglass dispersion |16) . 

For the used parameters as the frequency exceeds 
10 meV, the commensurate maximum splits into incom- 
mensurate peaks [Figs.[2jb)-(e)]. In the range 15 meV < 
oj < 70 meV peak intensities on the diagonals of the 
Brillouin zone exceed somewhat their value on the zone 
edge. For larger frequencies the susceptibility becomes 
nearly isotropic around Q. As seen from Fig. [21 for 
moderate frequencies x"(kw) consists of several closely 
spaced peaks. Apparently, these peaks merge into one 
broad maximum due to crystal inhomogeneity, as ob- 
served in experiments of Refs. OEM]. 

The peaks in Fig.[2£b)-(e) correspond to spin excita- 
tions. Susceptibility (fl) coincides with the spin Green's 
function, and the vanishing real part of its denominator 

D r (koj) = oj 2 - cjRcTT(kw) - cj£ = (12) 

at a small imaginary part Di(koj) = wlmiT(kw) defines 
the dispersion of these excitations. Figure [4] demon- 
strates a graphical solution of Eq. (TT21 for a wave vec- 
tor near Q. The real part has three zeros one of which 
falls into the region of large spin-excitation damping. 
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Fig. 4 The frequency dependencies of the real D r (the blue 
solid line) and imaginary Di (the red dashed line) parts of 
the denominator in Eq. JTJ for k = (0.97tt, 0.97tt), 5 = 0.001 
and A = 4.8 meV. 
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Fig. 5 The usual (black squares) and supplementary (red 
circles) dispersions of spin excitations. Parameters are the 
same as for Fig.[JJ Horizontal dotted lines indicate frequencies 
in panels of Fig. [2] The wave vector varies along the zone 
diagonal, k = (k,k). 



Therefore, it is not seen in x"(ku;). Two other zeros 
correspond to maxima in Fig. Etc)-(e). The zero with a 
lower frequency is similar to the respective zero of the 
real part of the susceptibility denominator in the hole- 
doped case [TB] - with some correction due to ReTT(kw) 
it is close to uj^. The zero at uj w 58 meV arises due to 
the region of anomalous dispersion in D r (kuj), which is 
related to the sharp minimum in Z)i(kw). Both pecu- 
liarities are consequences of the folded electron disper- 
sion in this class of crystals. Small equi-energy contours 
shown in Fig. [3] and dominant contribution of states on 
these contours into 77(kw) are inherent in these crystals 
for moderate frequencies. This leads to sharp minima 
in Dj(kw) when the wave vector k connects two such 
contours for a selected frequency. 
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Figure [S] demonstrates dispersions of two zeros of 
Eq. (fT2"]l . which fall into regions of small spin-excitation 
damping in Fig. |4l The locations of these zeros are close 
to the positions of maxima in Fig. EJc)-(e). As men- 
tioned above, the lower curve is close to the dispersion 
uj = cjk, and in this respect, it is similar to the spin- 
excitation dispersion in p-type cuprates. In the latter 
crystals, this dispersion forms the up-directed branch 
of the hourglass dispersion of susceptibility maxima. 
The difference between these compounds and the con- 
sidered crystal is in the lowest frequency of the spin- 
excitation dispersion ui r . For the parameters of Fig. [5] 
it is approximately equal to 6 meV, an order of magni- 
tude smaller than in moderately doped p-type cuprates. 
Besides the usual spin-excitation branch, in the consid- 
ered crystal there is a supplementary dispersion, which 
is connected with the second zero at higher frequencies. 
As mentioned above, its appearance is related to folded 
electron bands in n-type cuprates. 

Thus, there are some similarity in the dispersions of 
susceptibility maxima in the two groups of cuprates - 
one of the two up-directed branches in n-type crystals 
is similar to the upper hourglass branch in the hole- 
doped compounds. In the electron-doped cuprates the 
maxima of this branch are usually more intensive than 
peaks of the supplementary branch (see Fig. This 
latter branch is the first difference between the two 
groups. The second difference is the behavior of x"(kw) 
for uj < Lu r , where Eq. (lT2"j) has no solutions and, there- 
fore, there are no spin excitations [Fig.^a), the lowest 
horizontal line in Fig. Hy. As mentioned above, in this 
frequency range maxima in the momentum dependence 
of the susceptibility are caused by peaks in — ImTT(koj) 
in the numerator of the expression for x"(ko;). In p- 
type cuprates these peaks are at incommensurate mo- 
menta, which accounts for the down-directed branch of 
the hourglass dispersion |16j . In this branch, with de- 
creasing lo maxima of x"(ko>) move away from Q. In n- 
type crystals — ImTT(koj) peaks at the AF momentum 
due to electron bands folded across the AF Brillouin 
zone border. As a consequence in this case the suscep- 
tibility dispersion has the shape of a cone with the apex 
point at uj = 0, k = Q (see Fig. [JJ). 

The frequency dependence of the susceptibility at 
the AF wave vector is shown in Fig. [5J In panel (a) the 
calculated susceptibility is compared with experimental 
results from Ref. [TU]- In these calculations, the param- 
eter S was chosen to fit the location of the maximum to 
the experimental one. In the calculated susceptibility, 
the maximum approximately coincides with the bottom 
of the usual spin-excitation branch in Fig. [SJ The bot- 
tom of the supplementary dispersion is seen as a weaker 
maximum at ui « 45 meV. Although the measurements 




to (meV) 

Fig. 6 The frequency dependence of the susceptibility at 
k=Q. (a) Experimental data in PLCCO with x = 0.12 and 
T c = 24 K at T = 2 K \W\ (circles) and the calculated x" for 
5 — 4 x 10~ 3 and A = 9.6 meV (the curve), (b) the calculated 
susceptibility for <5 = 5 x 10 -4 and A = 4.8 meV. 



[10] were limited to ui < 30 meV, some growth in the 
experimental susceptibility may be also observed near 
the upper edge of this frequency window. In the local 
magnetic susceptibility measured in Ref. 4^ a shoulder 
or a maximum may be revealed at ui ~ 30 — 45 meV, 
which can be also related to the supplementary spin 
excitations. In recent work [TT] a PLCCO crystal with 
x = 0.12 and T c — 21 K was investigated. In com- 
parison with the T c — 24 K sample of Ref. [TU] the 
T c = 21 K sample has smaller superconducting gap 
A « 5.5 meV and is presumably located much closer 
to the boundary of the long-range AF order. Therefore, 
this crystal has to be characterized by a much smaller 
parameter S than the T c = 24 K sample. Two maxima 
at 2 and 9.5 meV were observed in the local suscepti- 
bility at T — 2 K. Calculations in Fig. [51b) carried out 
with smaller than in panel (a) values of the parame- 
ters A and 6 reproduce qualitatively this experimental 
spectrum. Maxima in this figure correspond to the two 
branches of the spin-excitation spectrum. Notice that in 
Ref. |llj a maximum at ui = 2 meV was also observed 
in the T c = 24 K sample. However, this maximum was 
much weaker than the peak at 10.5 meV and the max- 
imum at 2 meV in the T c — 21 K sample. I could not 
obtain such kind of spectrum. It is conceivable that the 
weak maximum is connected with sample inhomogene- 
ity 
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4 Conclusion 

In this work, the two-dimensional t-J model and the 
Mori projection operator technique were used for inter- 
preting some peculiarities of the magnetic susceptibil- 
ity in the optimally doped PLCCO crystal. The electron 
dispersion derived from photoemission data was applied 
in the calculations. The dispersion is characterized by 
the band folding across the antiferromagnetic Brillouin 
zone border. This band folding plays a central role in 
the formation of the low-frequency commensurate re- 
sponse and in the appearance of the supplementary 
branch of spin excitations. The origin of these pecu- 
liarities is related to small equi-energy contours of the 
electron dispersion, which are the consequence of the 
band folding. For small frequencies, the nesting vec- 
tor of the contours is the antiferromagnetic momentum 
Q = (w, ir). Therefore, for such frequencies the imagi- 
nary part of the polarization operator — ImTT(kcj) peaks 
sharply at Q. This operator is in the numerator of the 
formula for the susceptibility x"(kw), and at a weaker 
momentum dependence of the denominator, it controls 
the behavior of the susceptibility, leading to the com- 
mensurate response. For moderate frequencies, the peak 
in — ImTT(kcj) splits and shifts to incommensurate wave 
vectors. An anomalous dispersion in the real part of 
the polarization operator, which is connected with the 
peak, leads to the appearance of three poles in the spin 
Green's function. Two of these poles fall into regions of 
a small spin-excitation damping and correspond to the 
usual and supplementary spin-excitation branches. The 
latter branch nested into the former leads to a comb of 
closely spaced peaks in momentum cuts. Presumably 
this structure is not resolved in experiments, being seen 
as a broad commensurate peak up to ~ 100 meV. It is 
very likely that the supplementary branch was already 
observed in the recent measurements of the frequency 
dependence of the susceptibility. 

Comparing with moderately doped p-type cuprates 
a common element in dispersions of susceptibility max- 
ima can be revealed - it is branches of usual spin exci- 
tations. However, in the PLCCO there exists also the 
supplementary branch nested into the usual one. Be- 
sides, in the PLCCO the bottom frequency of the usual 
branch uj r is an order of magnitude smaller than in 
hole-doped crystals. For cu < u> r the behavior of the 
susceptibility in the two groups of compounds is drasti- 
cally different. In both groups this behavior is governed 
by — ImiT(kcj) in the numerator of the susceptibility 
formula. However, if in the PLCCO this quantity peaks 
at Q, which leads to the commensurate low- frequency 
response, in hole-doped crystals — ImiT(kw) reaches a 



maximum value at incommensurate momenta, and the 
low-frequency response is incommensurate. 
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